Asymptotic Analysis of the Paradox 
in Log-Stretch Dip Moveout 
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Abstract 

There exists a paradox in dip moveout (DM0) in 
seismic data processing. The paradox is why Notfors 
and Godfrey's approximate time log-stretched DM0 
can produce better impulse responses than the 
fuU log DM0, and why Hale's f-k DM0 is correct 
although it was based on two inaccurate assumptions 
for the midpoint repositioning and the DM0 time 
relationship? Based on the asymptotic analysis 
of the DM0 algorithms, we find that any form of 
correctly formulated DM0 must handle both space 
and time coordinates properly in order to deal with 
all dips accurately. The surprising improvement 
of Notfors and Godfrey's log DM0 on Bale and 
Jakubowicz's full log DM0 was due to the equivalent 
midpoint repositioning by transforming the time- 
related phase shift to the space-related phase shift. 
The explanation of why Hale's f-k DM0 is correct 
although it was based on two inaccurate assumptions 
is that the two approximations exactly cancel each 
other in the f-k domain to give the correct final result. 

Citation detail: X. S. Yang and B. Zhou, Asymp- 
totic analysis of the paradox in log-stretch dip move- 
out, Geophys. Res. Lett., 27(3), 441-444 (2000). 



INTRODUCTION 

Dip Moveout (DM0) technique has been widely used 
in seismic data processing over the past decade and 
many different algorithms have been developed [Hale, 
1984; Bale and Jakubowicz, 1987; Notfors and God- 
frey, 1987; Liner, 1990; Gardner, 1991; Black et 
al, 1993; Zhou et al, 1996]. Hale's DM0 is ac- 
curate for all reflector dips and has become an in- 
dustrial standard, but it is computational expensive 
and temporally non-stationary. In all the methods, 
a logarithmic time-stretching technique [Bolondi et 



al., 1982] is widely used, probably due to its com- 
putational efficiency and easy implementation. The 
impulse response produced by Bale and Jakubow- 
icz's full log-stretch operator (hereafter referred to as 
Bale's full log DM0) in the frequency-wavenumber (f- 
k) domain is inaccurate although no approximations 
were made in the mathematical formulations. The 
impulse response is surprisingly improved by Notfors 
and Godfrey's approximate log-stretch scheme (here- 
after referred to as Notfors' log DM0). Bale's inaccu- 
rate DM0 impulse responses reveal that Hale's DM0 
derivation had some inappropriate assumptions. To 
improve the DM0 impulse responses, an exact log- 
stretch DM0 was derived by transforming Hale's 
time log-stretch impulse response into the Fourier do- 
main [Liner, 1990] . Black et al. [1993] provided 
another derivation of the more reliable f-k DM0 by 
treating the reflection-point smear from dipping re- 
flectors correctly. Based on Black's DM0 relation- 
ship, Zhou et al. [1995, 1996] presented a log-stretch 
f-k DM0 which is similar to Gardner's DM0 [Gard- 
ner, 1991] but from a derivation with more direct 
physical insight into the DM0 process. All these 
DM0 schemes were developed to try to generate bet- 
ter impulse responses and to improve the computa- 
tional speed, but have different corrections of Hale's 
subtle flaw and improving Bale's inaccurate impulse 
responses. 

Despite extensive studies and the routine use of the 
DM0 algorithms, there still exists a paradox which 
has not yet been explained satisfactorily. The para- 
dox is why Notfors' approximate log DM0 can pro- 
duce better impulse responses than Bale's full log 
DM0, and why Hale's f-k DM0 is correct although 
it was based on two inaccurate assumptions. Liner 
pointed out [Liner, 1990] that Bale and Jakubow- 
icz's DM0 derivation implicitly assumes that the 
Fourier transform frequency in the log-stretch domain 
is time- independent. Black et al. attributed Hale's 
subtle flaw to the lack of midpoint change and the 
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consequently improper treatment of the reflection- 
point smear from dipping reflectors. Zhou et al. 
suggested that the approximations from the non- 
repositioned midpoint and from the incorrect time 
relationship before and after DMO may counteract 
each other in the f-k domain to make the final result 
correct, but they may not cancel each other in the log- 
stretch f-k domain. In this paper, we analyse these 
viewpoints and compare DMO algorithms to find out 
what is really responsible for the inaccurate impulse 
responses of Bale's full log DMO scheme. For con- 
venience in the following discussion, only common- 
offset DMO forms will be analysed although simi- 
lar methodology can be also applied to common-shot 
DMO analysis. 

ASYMPOTIC ANALYSIS OF 
PHASE SHIFT FUNCTIONS 

We assume a constant-velocity medium and follow 
the notations used earlier [Hale, 1984; Notfors and 
Godfrey, 1987; Black et al, 1993; Zhou et al, 1996]. 
The heuristic DMO mapping can be expressed by 

PQ(to,Xo,h) ^ Pn(tn,Xn,h), (1) 

where Pn{tm Xn, h) is a normal- moveout (NMO) cor- 
rected input section and PQ{tQ,XQ,h) is a zero-offset 
DMO output section, h is the half offset, t„ and a;„ 
are time and space in the NMO-corrected domain re- 
spectively, and tf) and xq are the time and space in 
DMO-corrected domain, respectively. The time log- 
stretch transform pair [Bolondi et al, 1982; Bale & 
Jakubowicz 1987] 

T = \n{t/tc) and t = t^e^ , (2) 

transforms the time coordinates from {to,tn) to the 
log-stretch variables (to,t„). tc is the minimum cut- 
off time introduced to avoid the singularity of the 
logarithm at zero. 

The DMO processing algorithm can be divided into 
the following five steps [Liner, 1990; Gardner, 1991; 
Hale, 1991]: 1) transform NMO corrected input data 
Pn(tn,Xn,h) to Pn{Tn,Xn,h) by the time log stretch 
relation ([2]), 2) 2-D FFT of P„(t„, x„, h) over t„ and 
Xn to Pn{^, k, h), 3) use the phase factor for filtering 
in the (57, k) domain to get Po(^^: k, h), 4) inverse 2-D 
FFT of PQ{n,k,h) over fl and k to Po{To,XQ,h), and 
5) inverse log stretch of Po(''o, a^O: ^) by the inverse 
relation of equation ^ to get the output Po{tQ,xo, h). 
Of all the available log stretch DMO algorithms, the 
same transform relation © is used in steps 1) and 5), 
and the two-dimensional FFT transform pair is also 



used in steps 2) and 4). The main difference between 
algorithms is thus the multiplying phase factor in step 
3). The following discussion will therefore focus on 
the analysis of the geometry-determining phase factor 
in different log-stretch DMO algorithms. As we will 
see, all the phase shift functions depend mainly on 
the variable 



Mathematically speaking, ^ can vary from (zero dip 
at the centre of the DMO ellipse) to co (very shallow, 
steeply dipping). 

Bale's Full Log DMO 

From Hale's DMO formulation. Bale and Jakubowicz 
[1987] derived the full log stretch DMO 

Po(f^,fc,/i) =e-*5™n(i-C^)p„(r!,fc,/i), (4) 

However, this algorithm does not produce good im- 
pulse responses. Liner [1990] attributed this feature 
to the fact that Bale and Jakubowicz's derivation im- 
plicitly assumes that f2 is not time-dependent. In 
fact, Bale's phase shift function 

^F = -\^Hi-e), (5) 

is only physically meaningful when ^ < 1 (near zero 
dips). It becomes complex if ^ > 1 (shallow, steep 
dips) which is responsible for the inverted Gaussian 
shape. 

Notfors' Log DMO 

An approximation to ([4]) was derived by Notfors and 
Godfrey [1987] 

Po(f^, fc, h) = e'"(^i+«^-i)p„(f], fc, h), (6) 

which is based on two approximations: (1) the inde- 
pendence of and T„, and (2) (hk/VlY << 1. With 
these approximations. Bale's phase shift function be- 
comes Notfors' phase shift function 

$Ar = r!(v/TTe2-i), (7) 

which is meaningful for any value of ^ (all reflector 
dips). This derivation implicitly extends the defini- 
tion range of ^ although the approximations are based 
on ^ << 1. It is this second approximation that 
makes Notfors' log DMO produce better impulse re- 
sponses than Bale's full log DMO @. 
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Liner's Log DMO 



Case 1: ^ << 1 



The exact log stretch DMO derived by Liner [1990] For C ^ 0, equations ©, ©, ^ aU become 
is 



PQin,k,h) = 



Pn{n,k,h), 



where 



Ps = Vs/h, 

A. - iln(l - / 
and the stationary point is 



h 



(8) 



(11) 



F.N,E 



(14) 



which means that all these log-stretch DMO algo- 
rithms are asymptotically equivalent near zero dips 

(9) 

— > 0) and will yield virtually the same impulse re- 
sponse geometry near the center of the DMO ellipse. 

(10) 

Case 2: ^» 1 

As the reflector dips become shallow and steep (for 
^ >> 1), expression ^ will be invalid, but (O and 
(|13p can be asymptotically approximated as 



Liner's derivation introduces log stretch variables 
into Hale's {t, x) elliptical response rather than into 
the NMO equation, as was done earlier by Bale 
and Jakubowicz [1987]. The formula ([S]) does yield 
the correct impulse response geometry [Liner, 1990; 
Zhou et al, 1996]. 
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Zhou's Log DMO 



Zhou et al. [1996] presented an accurate log stretch 
DMO 



and 



.m-i|), 



(15) 



(16) 



respectively. Using ^ — )• as ^ — ?► oo, we can write 
([15]) and ([TB|) as a single asymptotic form 



N,E 



(17) 



Po{n,k,h)^e' 



(12) 

which is equivalent to Gardner's results except for a 
sign difference on the phase term due to a different 
definition of the 2-D Fourier transform in Gardner's 
formulation. Zhou et al. [1996] derived their result 
based on the DMO relationships given by Black et 
al. [1993]. It is easy to check that the phase shift 
in equation is identical to that in equation ([5]). 
Therefore, Zhou's log DMO and Liner's log DMO will 
produce the same impulse response geometry. The 
only difference is that Zhou's DMO will yield slightly 
larger amplitudes than Liner's DMO because ^ has 
an amplitude factor 1/(1 -I- Pl)^/'^ < 1. For con- 
venience in the following analysis, we will refer to 
Liner's log DMO and Zhou's log DMO as the exact 
log DMO, and will use the simpler phase shift expres- 
sion 

<fE = ln\^iTW-i-HliVi + ^+i)l (13) 

which is accurate for all reflector dips. 

All the phase shift functions ([S]), (O and are 
monotonically increasing, therefore, we discuss only 
two asymptotic cases. 



which means that Notfors' log-stretch DMO and the 
exact log-stretch are asymptotically equivalent but 
they are not approximately equal due to the term 
'O(^) in (HH). This implies that the exact log 
DMO is able to deal with all reflector dips correctly, 
but Notfors' log DMO mishandles the shallow, steep 
dips although it is an improvement on Bale's DMO 
scheme. 



THE EXPLANATION OF THE 
PARADOX 

Based on the above analysis, we can explain the para- 
dox of why Notfors approximate log DMO improved 
Bale's full log DMO, and find out the real flaw in 
Hale's formulation. The f-k integral for common- 
offset in Hale's DMO [1984] is 

Po{uj,k,h)= [ f Po{to,x„,h)e'^^'°-''^»Ux,,dtn 



where 



A-fcx„)p^ (t„ , Xr, , /l)dx„d<„ , (18) 



^-A/i + (r#^)^ = AA + (/^)^- 

tn dXQ V tnt^ 



(19) 



3 



and the DM0 mapping ^ and the following relations 
of time and space coordinate are used 



to = Atn and 



xo 



(20) 



which are not strictly correct and make the DM0 op- 
erator amplitude-unpreserved. The similar time and 
space relationships were used by Bale and Jakubow- 
icz [1987] in their full log DM0 derivation. Black et 
al. [1993] derived a new amplitude preserving DM0 
integral 



Po(w,fc, /i) 



2^2 - 1 

Ti ' 



Po{to,xo,h)e'<^^*'>-k.o)d^^dt^^ 



(21) 

by substituting ([1]) and basing on Black's correct re- 
lationship of time and space coordinates before and 
after DM0 processing 

to ^ -r and xo=Xn~-r- , (22) 

A A tn (JJ 



where the definition of A is the same as in ([191) . 
From equation (|22p . we can easily derive the DM0 
impulse ellipse (Zhou et al., 1995) while the rela- 
tionships in equation ((20)) give an impulse point. As 
A > 1, Black's correct relationship (P^ implies the 
inequality 

to < tn. (23) 

It is clearly seen that phase shifts in Black's DM0 
(PT|) and in Hale's DM0 ^ are identical [Black 
et al., 1993] due to the identical phase factor 
exp[i(ujtnA—kxn)]- This implies that the two approx- 
imations resulting from the inaccurate relationships 
(1^(71) exactly cancel each other in the f-k domain, and 
thus make the final result correct. However, the coun- 
teraction of the two approximations does not give ex- 
act cancellation in Bale's log-stretch DM0 algorithm 
because only the time coordinate is stretched and the 
space coordinate remains unchanged. Therefore, we 
expect that the phase shift in Bale's log DM0 is es- 
sentially due to the log-stretch traveltime relationship 



1. 



To 



-Mi-e). 



This implies the incorrect inequalities 
To > r„ and > tn 



(24) 



(25) 



which contradict the correct relation (1^^ . Notfors 
and Godfrey [1987] used the approximate version of 



which yields the same incorrect inequality as in (j25l) . 

From the general forms of log-stretch DM0 algo- 
rithms, we know that the phase shift function 



$ = k{Xn - Xo) - VL[Tn - Tq), 



(27) 



consists of two parts: a space-related phase shift 
k[xn — Xo) and a time-related phase shift 0.{Tn — tq). 
The asymptotic analysis in the previous section re- 
veals that the space-related phase shift is always 
greater than the time-related phase shift. Compar- 
ing the phase shift functions (O and ([T^ , we see that 
Notfors approximate log DM0 virtually transformed 
the time-related phase shift into its space-related 
counterpart, and thus equivalently repositioned the 
midpoint 



Xo 



xn + j{i-^i+e). 



(28) 



To - T, 



(26) 



where the second term on the right hand side is simi- 
lar to the expression ()11|) of the stationary point in 
Liner's exact log DM0. It is this equivalent midpoint 
repositioning that makes Notfors' approximate log 
DM0 overcome Hale's subtle flaw and subsequently 
produces a better match to the geometry of the im- 
pulse response DM0 ellipse, but the Notfors' wider 
geometry of the DM0 ellipse shows that the shallow 
steep dips are still improperly treated. 



CONCLUSIONS 

The analysis of different log-stretch DM0 algorithms 
has shown that any form of correctly formulated 
DM0 must handle both space and time coordinates 
properly in order to deal with all dips accurately. 
The log-stretch DM0 algorithms presented by Liner 
[1990] and Zhou et al. [1996] are able to treat all 
reflector dips correctly. Bale and Jakubowicz's full 
log DM0 is only approximately correct for the near 
zero dips, and it becomes singular for shallow, steep 
dips. Notfors and Godfrey's approximate log DM0 
is a great improvement in handling all dips, but its 
impulse response still departs from the exact DM0 el- 
lipse, especially for shallow steep dips. The surprising 
improvement of Notfors and Godfrey's log DM0 on 
Bale and Jakubowicz's full log DM0 was mainly due 
to the equivalent midpoint repositioning by equiv- 
alently transforming the time-related phase shift to 
the space-related phase shift. The explanation of why 
Hale's f-k DM0 is correct although it was based on 
two inaccurate assumptions is that the two approxi- 
mations exactly cancel each other in the f-k domain 
to give the correct final result. 
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